NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



DOE/NASA/2674-79/8 
NASA TM-7930S 


EFFECT OF VELOCITY OVERSHOOT 
ON THE PERFORMANCE OF 
MAGNETOHYDRODYNAMIC 
SUBSONIC DIFFUSERS 


Mahesh S. Greywall 
Wichita State University 

and 

J, Marlin Smith 

National Aeronautics and Space Administration 
Lewis Research Center 


Work performed for 

U.S, DEPARTMENT OF ENERGY 
Energy Technology 
Magnetohydrodynamics Division 

(NASR-TH-79305) EFFECT OF VELOCITY 
OVERSHOOT ON THE PERFORMANCE OF 
MAGNETOHYDRODYNAMIC SUBSONIC DIFFUSERS 
(NASA) 10 p HC A02/MF A01 CSCL 201 Dnclas 

G3/75 il6445 

Prepared for 

Eighteenth Aerospace Sciences Meeting 

sponsored by the American Institute of Aeronautics and Astronautics 
Pasadena, California, January 14-16, 1980 



N80- 14922 


NOTICE 


This report was prepared to document work sponsored by 
the United States Government. Neither the United States 
nor its agent/ the United States Department of Energy, 
nor any Federal employees, nor any of their contractors, 
subcontractors or their employees, makes any warranty, 
express or implied, or assumes any legal liability or 
responsibility for the accuracy, completeness, or useful- 
ness of any information, apparatus, product or process 
disclosed, or represents that its use would not infringe 
privately owned rights. 


DOE/NASA/2064-79/8 
NASA T/Vl-79305 


EFFECT OF VELOCITY OVERSHOOT 
ON THE PERFORMANCE OF 
MAGNETOHYDRGDYNAMIC 
SUBSONIC DIFFUSERS 


Mahesh S. Greywall 
Wichita State University 
Wichita, Kansas 67208 

and 

J. Marlin Smith 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Prepared for 

U.S. DEPARTMENT OF ENERGY 

Energy Technology 

Magnetohydrodynamics Division 

Washington, D. C. 20545 

Under Interagency Agreement EF -77 -A -01 -267 4 


Eighteenth Aerospace Sciences Meeting 

sponsored by the American Institute of Aeronautics and Astronautics 
Pasadena, California, January 14-16, 1980 


THE EFFECT OF VELOCITY OVERSHOOT OH THE PERFORMANCE OF 
MAGHETOHYDRODYHAMIC SUBSONIC OIFFUSERS 

Mahesh S. Groyv/all 
Wichita State University 
Wichita, Kansas 67208 

and oia IClltj 

J. Marlin Smith " \ . - i To VOOB- 

NASA Lov(is ReS(?arch Center . . 

Cleveland, Ohio 44135 ■ 


Abstract 

The effect of velocity overshoot, a 
velocity distribution encountered in 
magnetohydrodynamic channels with peak velocity 
near the wall instead of at the center, on the 
performance of two-dimensional subsonic 
diffusers was investigated, For the simplified 
conditions analysed in this paper it was found 
that the ratio of peak-to-centerline velocity 
increases along the diffuser for a diffuser half 
angle greater than some critical value, It was 
also found that irrespective of the accompanying 
inlet temperature distribution the wall shear 
stress and the v/all heat flux is substantially 
larger when the inlet velocity profile has an 
overshoot than that for a fully 
developed turbulent inlet velocity profile. 

I . Introduction 

In ordinary channel flow the fluid is acted 
upon by a pressure force which is essentially 
constant across the channel and a viscous force 
which is zero at the center and maximum at the 
walls. Governed by these forces the velocity 
profile approaches a shape with maximum at the 
center (parabolic for laminar flows and flatter 
profiles for turbulent flows). In 
magnetohydrodynamic (MHD) power generator flows, 
the flow is also acted upon by a'j,xB Loreiitz 
force. This is a retarding force^wTien power is 
being extracted from the channel. Along the 
electrode wall, the J.x^ force is the same as in 
the free stream due ^"the requirement of 
current conservation. However, along the 
insulating wall the rapid decrease in electrical 
conductivity due to the large temperature 
differential between the free stream and the 
wall results in a similar decrease in the j,xB 
force. As a result, the flow is slowed down” 
more rapidly in the free stream than in the 
insulator wall boundary layer, which, depending 
upon flow conditions and the magnitude of the 
JxJI force, can cause the maximum velocity of the 
TTow profile to shift from the centerline into 
the boundary layer. This occurrence of a 
velocity maximum near the wall is termed 
velocity overshoot. These velocity overshoots 
have been theoretically predicted (references 1 
and 2) and experimentally observed (ref. 3). 

In the MHD power train a diffuser is 
attached to the exit of the generator cliannel. 
This paper theoretically investigates the effect 
of the velocity overshoot on the performance of 
the subsonic diffuser. The study is divided 


into two parts. In part one, given in section 
II, the evolution of an overshoot velocity 
distribution is Studied in a plane 
two-dimensional diffuser as a function of 
diffuser divergence angle. In part two, given 
in section III, the diffuser performance for 
velocitv overshoot is compared to that for fully 
developed inlet velocity profile. 

The computations were carried out using a 
newly developed method iref. 4) to calculate 
two-dimensional momentum and energy transport 
for a viscous compressible flow. In this method 
the flow in the duct is partitioned into finite 
streams. The difference equations are then 
obtained by applying momentum and energy 
conservation principles directly to the 
individual streams. For simplicity in our 
initial analysis, all calculations have been 
carried out for air assuming a perfect gas with 
constant specific heat and an inlet static 
pressure of 1 atm. 

Evolution of Velocity Ove rshoot as^^_a_ 

Function df'Dfff user A rTgT e 

In initial calculations utiliz'ng a 
two-dimensional diffuser with a u.-o oeg/ee half 
angle, an unexpected phencinena was observed. It 
was observed that the ratio of the peak velocity 
(velocity of the overshoot) to the centerline 
velocity increased with distance down the 
diffuser. This v;as unexpected since it was 
initially thought that viscous forces would tend 
to dampen the affect of the velocity overshoot. 

In order to more closely examine the above 
phenomena a study was carried out to investigate 
the effect of friction on the damping of the 
overshoot at various diffuser half angles. 
Furthermore in order to minimize any influence 
that heat transfer might have upon the 
phenomena, the calculations were carried out 
with the inlet centerline and wall temperatures 
set equal to 200DK. The inlet Mach number based 
on the centerline velocity was set equal to .8 
which is in the neighborhood of the expected 
exiting Mach number from a subsonic MHD power 
generator. The inlet velocity profile (see 
figure 1) was choosen as a convenient 
mathematical simulation of the velocity 
overshoot profile calculated in reference 1. 

Figure 1 shows the evolution of the 
overshoot velocity distribution in a 
two-dimensional diffuser with a two degree half 
angle. In this and the other figures the 


following nomonclatupo is usod: 

X, distance along the diffuser centerline, 

2 , distance from the diffuser wall, 

Hq, diffuser height at the inlet, 

tl, diffuser height at x. 

As previously discussed, figure 1 shows 
that while the velocity as a whole decreaseci 
down the diffuser, that the ratio of the 
peak-tO“centerlino velocity does not decrease, 

In fact this ratio increases with distance along 
the diffuser. The fiTspected reason for this 
somewhat unexpected behavior is that while 
viscous forces tend to smear out the peak, the 
rising pr essure in the decelerating flow 
TcTenTuaTes the velocity differences normal to 
the flow, i.G. a fast moving stream has to slow 
down less than a slow moving stream to overcome 
the same pressure rise. 

In order to check the validity of the above 
observation several cases were run with 
different diffuser divergence angles since if 
the above explanation of the phenomena is 
correct then at a small enough divergence angle 
the friction forces must eventually dominate and 
result in a damping of the overshoot. Figure ? 
shows the evolution of the overshoot inlet 
velocity distribution for zero divergence, angle, 
i.e. for flow through a constant area duct. We 
note that for this case the ratio of 
peak-to-centerline .velocity decreases with 
distance down the duct. In this case pressure 
falls along the duct and the decreasing pressure 
helps the viscous forces to smear out the 
velocity overshoot. Figure 3 shows the 
variation of the peak-to-centorline velocity 
ratio along the diffuser for various divergence, 
angles. In all cases the flow conditions at the 
inlet are identical. The leveling off of the 
growth in peak-to-centerline velocity ratio for 
two and three degree half angles is due to the 
leveling off of the pressure rise as the flow 
travels down the diffuser. 

m • Comparison of Diffuser Performance 
for V eTociFy^vershdot an?TuTT.y ' 
PeveTbliegTnTet Prdfi 

In actuality both the velocity and 
temperature profiles are specified by the MHD 
generator exit conditions. The degree of 
velocity overshoot, if any, and the effect upon 
the temperature profile are strongly dependent 
upon the operating characterisitics of the MHD 
generator. Therefore, to investigate the 
general character of the effects associated 
with velocity overshoot, three d'ifference 
representative combinations of jnlet velocity 
and temperature profiles were eximi"ned. To 
provide a common basis on which to compare the 
three different cases, the representative inlet 
conditions were selected for case 1 and 
conditions for the two other cases (cases 2 and 
3) were selected to match the mass flow rate and 
stagnation enthalpy of case 1. In all cases the 
inlet static pressure and wall temperature were 
1 atm and 1500K respectively. 

Case 1: Fully developed Velocity - Fully 
Developed Temperature - In this case the 


velocity u and temperature T at the inlet are 
calculated from the formula* 

u T - T„ z 1''^ 

D r.’ 

’'"c " ^w ^c 

where subscription e aiml w refer to the wall and 
centerline conditions rospectiyely. 

Typical numerical values for the exit of a 
MHD generator were selected for tho centerline 
velocity and temperature of 717.. F M/S and 200K 
respectively. Those correspond to a Mach number 
of 0.8. 

C ase P! Overslino t Velo ciFy - Fully , DAVAIQ^'I'^ 
teinpef a t~u r~o " raclT1¥g‘^ecTf^¥'HafrTr^d”the 
expli eft form of the overshoot velocity profile, 
a mathematically convenient approximation of the 
profile given in reference 1 was used. In this 
approximation the variation of velocity from the 
wall to the peak is taken as V'' and the peak 
velocity is fifty percent greater tlian the 
centerline velocity. Of even greater 
uncertainty is tho associated temperature 
profile. It is therefore the choice of tho 
profile that separates cases two and three, In 
this case the limit of a fully developed 
temperature profile was assumed as in the 
previous case, 

T - Tv, ^ z 1/7 

’’’c ■ fw 7c 

The inlet centerline velocity and temperature 
were 56R.7 M/S and 2000. 7K respectively to match 
the mass flow rate and stagnation enthalpy of 
case 1. 

C ase 3: Overshoot V eloc ity - Overshoot 
Temper ature - InTlTisTase the velocity profile 
is fdentTcal to that of case two hut a liniiting 
case in which the temperature profile is 
identical to the velocity is considered so that 
the temperature is given by: 

T - T„ u 

1"c " fw 

and the equivalent inlet centerline velccity and 
temperature of 566.7 M/S and 1892K 
respectively. 

Calculated Result s 

The inlet velocity and temperature profiles 
for the above three cases are presented in 
figures da and 5a. 

The diffuser flow for the three specified cases 
was two degree half angle, one meter in height 
at the inlet, and ten meters long. Results of 
the calculations are shown in figures 4 through 
8. Figure 4b shows the velocity profiles at the 
diffuser e.xit for the three different cases. We 
note, that for the conditions considered here, 
that irrespective of the accompanying 
temperature profile the velocity overshoot 
persists up to the end of the diffuser. Figure 
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5b shows the temperature profiles at the 
diffuser exit. We note that the temperature 
overshoot for the OVERSHOOT VELGiITY-OVERSHOOT 
TEMPERATURE case has substantially disappeared. 
The exit temperature profile for the OVERSHOOT 
VELOCITY-FULLY DEVELOPED TEMPERATURE case shows 
a sudden change in the slope near the wall. The 
point of sudden change corresponds roughly to 
the point of peak velocity and is possibly due 
to the compression work associated with slowing 
down this velocity peak. 

Figure 6 shows the pressure rise along the 
diffuser for the three cases. The Insignificant 
difference between the different cases Indicates 
that the pressure rise is almost solely due to 
the area change and hence the pressure recovery 
should be nearly Identical for the three cases. 
However, from figure 7, which shows the shear 
stress at the wall, we note that irrespective of 
the accompanying temperature profile the 
velocity overshoot case has two and a half to 
four times the wall shear stress compared to the 
fully developed velocity profile case. Tnus, as 
noted in references 1 and 'i, the velocity 
overshoot condition should be less susceptible 
to separation than profiles of the fully 
developed form. This indicates that the 
electrode wall which does not possess velocity 
overshoot will probabl;/ separate before the 
insulator wall unless it is designed with 
smaller divergence angle. 

Figure 8 shows heat flux to the wall for 
the three cases. The substantially larger heat 
flux to the wall in the case of OVERSHOOT 
VELOCITY-OVERSHOOT TEMPERATURE as compared to 
the FULLY DEVELOPED VELOCITY-FULLY DEVELOPED 
TEMPERATURE case is partly due to the higher 
temperature (temperature overshoot peak) in the 
vicinity of the wall and partly due to enhanced 
convection caused by the higiter velocity 
(velocity overshoot peak) near the wall. In the 
case of OVERSHOOT VELOCITY-FULLY DEVELOPED 
TEMPERATURE the increase in heat flux is 
primarily due to the increased convection. As 
we proceed down the diffuser the temperature 
overshoot, as seen from figure 5, gradually 
disappears and the heat flux to the wall for the 
two cases OVERSHOOT VELOCITY-OVERSHOOT 
TEMPERATURE and OVERSHOOT VELOCITY-FULLY 
DEVELOPED TEMPERATURE should approach each 
other. This is illustrated in the results shown 
in figure 8. The large heat flux to the wall 
exhibited by the velocity overshoot cases 
indicates that the velocity profile exiting the 
MHD generator will also have a serious effect 
upon the overall heat transfer design of the 
diffuser. 


IV. Concluding Remarks 

In this preliminary study of the effects 
upon diffuser performance resulting from the 
overshoot velocity profiles characteristic of 
the flow along the insulator walls at the exit 
of a MHD generator, a number of simplifying 
assumptions were made. A two-dimensional 
calculation was used, air was assumed to be the 
working fluid and in addition it was assumed to 
be a perfect gas, and a mathematically 
convenient representation of the velocity 


overshoot profile calculated in reference 1 was 
assumed. In actuality the working fluid will be 
combustion products and the degree of velocity 
overshoot, if any, will be determined by the 
performance characteristics of the MHD 
generator. Therefore the conclusions reached In 
this preliminary study are restricted by the 
limitations of the input conditions. A 
comparison between velocity overshoot and 
fully-developed profiles indicates the velocity 
overshoot condition to have significantly 
greater wall shear stress and hence should be 
less susceptible to flow separation than the 
fully-developed profile. The heat flux 
associated with the overshoot conditions Is 
however higher, and hence will have a serious 
effect upon the overall heat transfer design of 
the diffuser. An unsuspected result of this 
investigation was the amplification of the 
velocity overshoot for diffuser divergence 
angles greater than some critical value due to 
the dominance of pressure effects relative to 
the damping effects of viscosity. The effects 
of flow separation, half angle divergences, and 
total pressure recovery are presently being 
studied for the more realistic condition of a 
combustion gas working fluid and for a range of 
velocity overshoot magnitudes and profiles. 


REFERENCES 


1. Maxwell, C, D., Doss, E, D., Curry, B, P., 
Oliver, D. A., "Consideration of Three- 
Dimensional Effects in MHD Power Generators," 
Engineering Aspects of Magnetohydrodynamics , 
15th, University of Pennsylvania, 1976, 

pp. IX. 6.1 - IX. 6. 10. 

2. Rankin, R. R., and Eustis, R, H., "Boundary 
Layer Calculations for the MHD Insulator Wall, 
Engineering Aspects of Magnetohydrodynamics . 
15th, University of Pennsylvania, 1976, 

pp. 1.5.1 - 1.5.4. 

3. Rankin, R. R. , Self, S. A., and Eustis, R. H,, 
"A Study of the MHD Insulating Wall Boundary 
Layer," Engineering Aspects of Magnetohydro - 
dynamics . 16th, University of Pittsburgh, 1977 
pp. X.8,1 - X.8.11. 

4. Greywall, M. S., "Streamwise Computation of 
Duct Flows," Computer Methods in Applied 
Mechanics and Engineering (to be published). 


pv’T'"' ' 

OiUi: 


-{■■{ OF TllK 
. id POOR 


3 



OVERSHOOT VEL - OVKSHOOT TEft'iP. 








OVERSHOOT VEL - 
FULLY DEVELOPED TEMP. 
OVERSHOOT VEL - 



Figure 7. - Axial rail shear stress dlstribulion. ^ ^“butien. 



1 No 2 fco^mnitifti A««$ion No 3 Beop cnt't Ciiiios No 

NASA TM-79305 


0 ftspod Dji» 

f. Porlormm'j ‘Ajaniwiion &i3» 

0 Pcf(0fmm3 OfQ^nsfaiicfi Hcpn'«N'' 

E-257 _ ___ 

To Workl/nirNo 


4 T.i(« and Oobi.119 effect of velocity ovehshoot on 
THE PEHFORMANCE OF MAGNETOHYDKODYNAMie 
SUBSONIC DIFFUSERS 

TXiSr^lahesh S. GreYwaHTwIcFlt^^ 

Kansas, andJ. Marlin Smith, lewis Research Contor 


9 Performinj Ofsanuation Name end Addreu 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Olilo 4413D 


12 SponiorinB Ageney- Name and Addren 

U.S. Department of Eneray 
Magnetonydrodynamlcs Division 
Washington, D,C. 20545 


IB. Suppiemeniify Noies prepared under Interagency Agreement EF-77-A-01-2G74. Prepared for the 
Eighteenth Aerospace Sciences Meeting sponsored by the American Institute of Aeronautics and 
Astronautics, Pasadena, California, January 14-lC, 1980. 


16 Absifact 

The effect of velocity overshoot, a velocity distribution encountered in magnetohydrodynamie 
channels with peah velocity near the wall instead of at the center, on the performance of two- 
dimensional subsonic diffuserf wns bivestlgated, For the simplified conditions analyzed in this 
paper it was found tliat the ratio of poah-to-centerllne velocity increases along the diffuser for a 
diffuser half angle greater thajj some critical value. It was also found that irrespective of the 
accompanying bilet temperature distribution the wall slioar stress and tlie wall heat flux is sub- 
stantially larger when the inlet velocity profile has an overshoot than that for a fully developed 
turbulent inlet velocity profile. 


11 Contfict Of Gram No 


13 Typo Of Report aruj PcnoiJ Covered | 

Tocluiical Memorandum J 

14 Sponjonnp Aoency Cwlr Report NO.} 

DOE/NASA/2C74-79/8 f 


17. Key Wordi (Sugpi'sied by AuIhor(»ll 

Diffuser 

Magnetohydrodynamios 


IB Distribution Statomom 

Unclassified - unlimited 
STAR Category 75 
DOE Category UC 0OG 


19 Security Classif te( this report' 
Unclassified 


20 Security Classif. (ot this page! 

21. No of Pages 

Unclassified 




’ For sale by the National Technical Information Service, Spniielieltt, V'fBinia 22161 











